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FLUID  MECHANICS  AND  THE  TRANSPORT  PHENOMENA  1 

Many  of  the  current  research  problems  in  chemical  engineering  concern 
the  simultaneous  transfer  of  heat  and  mass  in  flow  system**  sometimes 
further  complicated  by  chonrwcal  reactions.  Such  flo-w  problems  axe  described 
by  the  !,equr.Uon*  ol  change*'  cf  fluid  mechanics.  These  basic  differential 
equations  form  the  starting  point  for  the  development  of  a number  of  topics 
of  direct  interest  to  chemical  engineering  science,  such  as;  (i)  the  solution 
of  laminar  problems  in  fluid  flow,  diffusion,  and  nc -*t  flow;  (ii)  the  study  of 
turbulent  flow  and  the  eddy  transport  processes;  (lii)  the  theory  of  boundary 
layers  and  films;  (iv)  the  study  of  flow  in  particulate  systems;  (v)  the  applica- 
tions and  limitations  of  analogies  between  mass,  momentum,  and  energy 
transport;  and  (vi)  the  dir-ensional  analysis  of  complex  flow  problems  wiih 
heat  and  mass  transfer,  which  cannot  be  solved  analytically.  Because  of  the 
paramount  importance  of  the  equations  of  change  in  basic  chemical  engineering 
studies,  we  summariee  these  equations  here  in  their  most  complete  form.  We 
further  point  out  the  major  relationchips  between  the  equations  of  change,  the 
flux  vectors,  the  transport  coefficients,  and  the  forces  between  molecules. 

The  equations  of  change  of  fluid  mechanics  comprise  the  equations  of 
continuity  for  each  of  the  chemical  species,  the  equation  of  motion,  and  the 


T Part  of  the  material  for  this  article  is  taken  from  "Molecular  Theory  of 
Gases  and  Liquids,''  by  Hir schfelder,  Curtiss,  and  Bird  (18).  In  this  recent 
book  (hereinafter  referred  to  as  MTGL)  the  equations  of  change  are  discussed 
in  detail,  and  their  applications  to  the  mathematical  theory  of  flam(?s,  shock 
waives,  sound  propagation,  and  detonations  are  presented. 
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equation  o i energy  balance.  From  the  first  of  these  equations  coarfA  the  basic 

' < M t '*  U;  i 

equations  for  diffusion,  from  the  second  comes  thr  basis  for  il»a  Study  of 
fluid  flow;  and  from  tba  last  equation  comes  the  starting  point  for  the  ftfudy 
of  "he  at  flow.  Most  of  the  chemical  engineering  unit  operations  involve  oim>  or 
more  of  these  processes.  Hence  advances  in  the  fundamental  undo?  standing 
of  the’ various  Irreversible  processes  occurring  in  the  c Ira  engineering 

opt  rations  depend  on  the  correct  interpretation  and  uee  of  the, equations  of 
change.  i ■ . , \ 

These  equations  of  change  are  coupled,  nonlinear  differential  equations, 
and  om  might- well  doubt  the  usefulness  of  such  a complex  set  of  equations 
beyond  the  pedagogical  value  of  providing  a cor-  act  formal  representation 
of  a number  of  related  fields.  Of  principal  interest  to  the  chemical,  engineer 
at**  the  solutions  fc»  these  equations.  These  solutions  can  be  divided  into 
three  categories;  analytical,  numerical,  and  experimental.  Analytical  ao,lu- 
tiona  can  be  obtained  for  very  simple  problems  only.  It  iff  ns  vox  the  lass  just 
these  simple  solutions  which  are  the  handy  formulae  for  every-day  calculations. 
Poise  uiile'a  law,  Stokes'  law,  Bernoulli's  equation,  and  the  formulae  for  mass 
transfer  coefficients  are  a few  examples.  Other  simple  solv&Uma,  which  are 
aomewkat  less  familiar,  describe  flow  near  .rotating. or  oacillMiwg.  discs,  flow 
of  plastics  through  tubes,  temperature  profiles  produced  by  viscous  heat  ef- 
fects, temperature  distributions  in  a hot  fluid  entering  a cold  pipe  fine  Gracia 
problem),  »r..d  the  thickness  of  falli-g  films  in  wettod-wall  tcwsXfl.  Numerical 
solutions  can  frequently  be  obtained  when  problems  ara  too  complex  for 
analytical  solutions.  The  current  developments  in  high  «p,;ed  computing 


techniquea  open  up  tremendous  new  frontiers  for  the  further  development  of 
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fluid  mechanics  am1  transport  processes  as  af.piaed  to  chemical  A3ngin.ee ring. 
Som.«  cxaicpk^  of  proolem#  attacked  by  numerical  method*  are  the  study -of 
flame  propagation,  the  formation  of  shock  waves,  flow  tha rough  mu  air's,  and 
flow  around  solid  objects.  Experimental  solutions  to  extremely,  complex  .. ,i . 
problems  form  the  basis  for  many  industrial  design  calculation*  in  diffusion^] 
operations,  heat  transfer,  and  chemical  kinetics.  The  experiment?.!  data  e&n 
usually  be  correlated  in  terms  of  several  dimensionless  groups,  The  minimum 
number  of  correlating  parameters  needed  can  be  deduced  directly  from  the 
equationsof  change  when  they  are  written  in  dimensionless  form.  Such  methods 
have  proved  quite  useful  in  analyzing  the  flow  around  cylinders,  the  evaporation 
of  droplets,  and  the  flow  patterns  of  convection.  In  fact,  the  dimensional  analy- 
sis of  the  basic  diffe rentlai  equations  gives  the  justification  for  the  friction 
factors,  drag  factors,  and  j-factors.  Hence  even  for  problems  with  turbulent 

? V-i 

flow  in  complex  geometries  the  equations  of  change  can  be  useful  even  though 
their  analytical  or  numerical  solution  would  be  impossible. 

In  this  paper  we  begin  by  giving  the  equations  of  change  in  thair  most 
general  form  in  terms  of  the  flux  vectors  Most  books  on  fluid  mechanics 


do  not  give  these  equations  for  fluid  mixtures  or  i'or  system*  with  temperature 

' ¥• 

gradients  Frequently  these  equations  axe  quoted  incorretHy  and  incompletely. 


it  is  hoped  that  the  summary  given  here  will  be  of  use.  IRTp  neat  proceed  to 


give  the  expressions  for  the  flux  vectors  in  terms  of  th*  transport  coefficients 


It  is  at  this  stage  where  approximations  are  generally  made  concerning  the 
mechanism  of  the  transport  processes  The  general  results  of  the  thermo- 
dynamics of  irre  ve  j sible  pro  eases  are  set  forth  here  Finally  we  discuss 
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the  transport  coefficient?  and  their  relation  to  inter  molecular  forces.  Since 
the  chemical  engineer  has  to  deal  with  fluids  of  widely  varying  physical 
properties,  this  aspect  of  fluid  mechanics  is  of  considerable  importance. 
Inasmuch  as  this  subject  has  already  been  treated  quite  extensively  in  a v w 
recent  review  article  (1),  only  the  main  points  of  the  .diecussioa  are  given 
here . 
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1 . The  Equations  of  Change  la  Terms  of  the  Fluxes 

The  behavior  of  a flowing  fluid  in  which  hast  and  vts.it  transfer  itt  taking 
place  is  described  by  the  equation  of  change,  along  with  the  thermal  and 
caloric  equations  of  state.  The  equations  of  change  are  differential  equations 
which  describe  the  changes  in  the  macroscopic  properties  of  tb£  fluid  (for 
example,  the  local  density,  stream  velocity,  and  temperature)  in  terms  of 
the  flux  of  ma«8,  the  flux  of  momentum,  the  flux  of  energy,  and  the  chemical 
kinetics.  Fir  st  we  discuss  these  equations  and  give  the  meaning  of  the  various 
terms  contained  therein;  then  rve  proceed  to  indicate  the  limits  of  their  ap- 
plicability . 

a.  Summary  of  the  Equations  of  Change 

Th»  basic  equations  of  continuity,  motion,  and  energy  balance  correspond 
to  the  fundamental  principles  of  conservation  of  mass,  momentum,  and  energy 
respectively.  They  have  been  derived  for  veiry  general  conditions  both  in 
classical  and  quantum  theory.  In  terms  of  the  fluxes,  the  at  equations  of 
change  for  a fluid  mixture  containing  nc  chemical  species  are: 

Equation  of  continuity  for  component  1 : 

Dfi/Dt.  = - (V-j,)  + K;  UU,™«v.  (I.U)f 

Equation  of  continuity  for  the  fluid  as  a whole: 

Df/Dt  = -p  lV-S)  (11b) 

t The  derivation,  of  the  equations  of  change  from  the  Boltzmann  equation  for 
dilute  monatomic  gases  is  given  in  MT G 1 , §7.2-c,d,  e.  The  analogous 
derivation  for  dense  monatomic  gases  le  given  in  MTGL,  $9.4d. 
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af  motion: 

^ A 

D*7Dt«  -0/i-XV-p)  + (i/f)2ft5  (12) 

“ '=  1 

&|Mtion  ©f  tcargy  balances 

»s 

EtJ/t*  - -0/fiXV-g)  - (1/pXp.Vir)  + (Vj.).Z  («)* 

In  these  equations  the  substantial  derivative  D/Ot  * A/A*  + (v.^JT)  represents 
the  time  rate  of  change  following  a fluid  element  which  is  Moving  with  a 
velocity  I . The  symbols  , p , and  q represent  the  flume  of  anese, 
raenwatm,  and  energy  with  respect  to  the  mass  average  velocity  _v  . These 
flams  may  be  broken  down  into  several  parte: 


ii  - p * f 

(1.4) 

£ - p£  ♦ f 

t 

(IS) 

That  Is,  tbs  mass  flu*  include*  components  due  to  gradients  in  the  canc  aUa  - 
tisa.  temperature,  and  pressure  and  an  additional  component  due  to  external 
forces.  The  momentum  flax  (pressure  tensor)  includes  s term  associated 
with  the  static  pressure  » sod  another  associated  with  viscous  stresses. 

» shpuld  be  mentioned  that  the  6 in  Kq,  (I.  i)  in  the  ttermodynanuc  interna! 
energy  (per  suit  massf,  and  does  not  include  the  kinetic  enargy  associated 
wills  too  stream  velocity  or  the  potential  energy  associated  with  the  external 
forces . 

jr  The  symbel  K*  ctands  for  the  net  rate  of  production  i*f  chemical  species  i 
bv  chomics)1  reaction,  in  units  of  gm/sec.. 
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An\l  ..he  energy  flux  Includes  terras  due  to  temperature  and  concentration 

i > • , . . 

gradients,  a term  which  accounts  for  the  energy  flux  accompanying  diffusion, 
and  the  energy  flux  due  tc  r adiation.  Explicit  expressions  for  these  various 
fluxes  axe  given  la  §2. 

According  to  Eq.  (1.  la)  the  mans  denBity  of  species  i contained  in  a 
fluid  element  changes  for  three  reasons,  the  first  term  on  the  right  indicates 
the  change  due  to  the  fluid  expansion;  the  second  term  accounts  for  the  change 
due  to  diffusion  processes  (this  includes  ordinary,  thermal,  pressure  and 
forced  diffusion);  and  the  third  term  represents  the  change  due  to  the  produc  - 
tion or  loss  of  species  i by  chemical  reactions.  'When  the  nc  equations 
given  in  Eq.  (1.  la)  are  added  together,  one  obtains  the  overall  equation  of 
continuity.  Eq.  (1.  lb).  The  latter  indicates  that  the  density  of  the  fluic  as 
a whole  changes  only  because  of  the  expansion  of  the  fluid  as  indicated  by  the 
single  term  on  the  right. 

According  to  Eq.  (1.2)  the  velocity  of  the  fluid  element  undergoes  a 

.i  .j  . * > ' .*•  • ' . . i»t* 

change  because  of  the  gradient  in  the  prsssur*  tensor  and  alzo  because  of 
the  external  forces  acting  on  the  various  chemical  species  praseut.  Since 
the  pressure  tensor  may  be  written  ae  the  sum  of  two  terms,  the  first  term 
on  the  right  hand  side  of  Eq.  (1.2)  contains  two  terms;  -(1/p)  Vp  which 
represents  the  acceleration  of  a fluid  element  because  of  a gradient  in  the 
static  pressure;  and  ~( 1 /p)(  V- p^)  which  gives  the  deceleration  of  the 
fluid  wLaaaont  because  of  internal,  friction  (viscous  streams). 

According  to  Eq,  (1.3)  the  Internal  energy  cha.igso  for  t-he  following 
reasons;  the  first  term  on  the  right  represents  the  change  due  to  energy 
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flux  (this  includes  thermal  conduction,  energy  flux  due  to  a concentration 
gradient,  energy  transport  by  molecular  diffusion,  and  energy  transport 
by  radiation);  the  second  term  Includes  both  the  energy  change  due  to  pV- 
work,  — ( 1 /jo)p(  V-  v)  , and  that  due  to  viscous  dissipation? , -(1  /pMp^ V^:  Vy)  ; 

and  the  third  term  describes  the  change  due  to  the  work  done  by  the  diffusing 
molecules  in  overcoming  tbe  external  forces. 

The  equations  of  charge  along  with  the  appropriate  boundary  conditions 
form  the  starting  point  for  any  mathematical  study  of  fluid  flow,  heat  transfer, 
and  diffusion.  For  the  complete  statement  of  the  problem,  one  also  needs  to 
specify: 

P " P fT.Pi , pz.,  • • • fnJ  (1.7) 

U - U (T,  p, , , . . . fin) 

which  are  the  thermal  and  c .loric  equations  of  state  respectively.  When 
explicit  expressions  are  introduced  for  the  flux  vectors  in  terms  of  the 
transport  coefficients,  then  there  appear  in  the  equations  of  change  the  co- 
efficients of  diffusion,  viscosity,  and  thermal  conductivity;  the  dependence 
of  these  quantities  on  the  temperature  and  the  mass -densities  of  tbe  various 
species  present  also  have  to  be  specified  in  the  solution  of  an  actual  problem. 

^ The  viscous  dissipation  terms  have  been  included  in  tbe  analytical  solution 
of  only  several  problems:  flow  in  capillary  tubes  by  Hausenblas  (17)  and 
by  Brinkman  (4);  flow  mconcentric  cylinders  by  Weltmann  and  Kuhn  (43) 
and  by  Blok  (2);  flow  near  a rotating  plate  by  Millaaps  aiul  Pohlhausen  (30), 
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The  equation  of  energy  balance  as  given  in  Eq.  ()..  3)  is  not  in  its  moat 
useful  (arm  tor  many  problems.  In  heat  transfer  problems,  for  example, 
it  is  considerably  more  convenient  to  rewrite  the  energy  equation  in  termu 
of  tho  temperature.  And  in  problems  on  sound  waves,  detonations,  and 
flow  through  nozzles,  it  is  more  appropriate  to  work  with  an  «qus.ii*>n  of 
change  for  entropy.  Furthermore  the  entropy  equation  is  important,  in 
the  formulation  of  the  thermodynamics  of  irreversible  processes. 

The  equation  of  change  for  temperature  may  be  obtained  by  using  F.q. 
(1.8)  to  rewrite  Eq.  (1.3).  When  no  assumptions  aj*e  iwqdb . ike  finhl 
result  is-. 

PC„ (DT/Dtj  - -(V-g)  - (j>t Var)  + (r-T(,Sp/&T))(V-i'0 

+ f.r(7-h)  + Kl][(Oi/m<)+(p-T(^^v;/K.i}]  (i-9) 

' n.  a 

+ L (jr  ?») 

in  whicfa  tlvs  derivative  (dp/dT)  is  taken  rt  constant  composition  «nd 
constant,  volume.  For  an  ideal  gas  the  terms  containing  (p  - T(dp/oT)) 
vanish,  and  Cor  a fluid  consisting  of  a single  chemical  fpooioa,  the  last 
two  terms  (these  involving  th«  summations  over  i ) need  not  be  considered. 

r » ‘ 

Hence  in  moat  heat  transfer  problems  with  fluid  flsr?,  one  hoc  tho  equa- 
tions cit  ceutiuruty  (Eq.  (i.  la)),  the  equation  of  motion  (Eq,  (1.2)),  and  the 
equation  of  eaergy  balance  in  terms  of  the  temperature  (Eq.  (I  .9)),  which 
have  to  be  solved  aktug  with  the  appropriate  initial  and  boundary  conditions. 
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Furthermore  one  needs  the  thermal  equation  of  state  (F.q.  (1. 7))  and  the 

A 

equation  for  Cv(T,jJi,  • - *pn  ) 'which  is  obtainable  from  the  caloric 
equation  of  state  in  Eq.  (1.8). 

The  equation  of  change  for  entropy  can  be  obtained  from  the  equations 
of  change  as  given  in  Eqs.  (1.1,  1.2,  and  1.  3)  and  the  differential  relation 
TdS  = dU  + pd(X/j>)  - (Uj^/mi)d{p^/jo)  , in  which  [ij  is  the  chemical 

W * 

potential5  (or  partial  molal  Gibbs  free  energy  5j).  The  equation  of  change 
for  entropy  can  be  written  in  the  form 

p(DS/Dt)  » — (V*ct)  + g <1.10) 


in  which  tr  is  the  flux  of  entropy  and  g is  the  rate  of  production  of  entropy. 


These  quantities  may,  after  considerable  minipulation,  be  written  in  the  form: 

r J ' >1.11/ 

— (1/T)  21  (jt- Ai)  - (i/T)ipW:  Vjr) 


O'  - 

3 ® 


(1/T)q(W  + X(Si/rr,i)j< 


(1-H)1 


t»l 


- (1/T)(c|(h).7  In  T)  +(|/T)2rj5 

S»> 


f All  partial  molal  quantities  used  in  this  article  are  on  a **por  molecule*1  basis. 
This  somewhat  unusual  convention  is  employed  because  of  the  fact  that  it  is 
customary  to  define  most  kinetic  theory  quantities  on  this  basis.  The  use  of 
the  rmr.dy  nasal  c quantities  per  molecule  then  eliminates  the  appearance  of 
Avogadro's  number  throughout  the  equations. 


I 


Here  the  total  energy  flux  has  been  arbitrarily  separated  into  two  parts. 


3 


w 


(ilia) 


thus: 


in  which 

q(T)  + q‘i>  + q(r» 


Sna 


«*# 


= Z ( Hi/m;)  ii 


(1.11b) 
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in  which  rB  and  Ya  are  the  net  rate  and  the  “chemical  affinity*  of  th^  e*1^' 
chssiical  reaction  (thebe  quantities  are  Jell  ied  in.  a footnote  below*  },  and 

4i  U:  n«.  ^ 

Ai*  (1/mi)Zo(^%*j)Tji»  Vpj  + (V^/mt)  Vf>  - ft  (Mi) 

The  quantity  A may  be  regarded  as  a driving  force  which  includes  the 
concentration  gradient,  the  pressure  gradient,  and  the  external  forces. 

Equation  (1.  11)  indicates  that  the  entropy  flux  is  made  up  of  two  components 
the  first  term  is  the  reversible  flux  of  entropy  due  to  heat  flow,  and  the  second 
term  is  the  flow  incidental  to  the  diffusion  processes.  Eq.  (1.  12)  indicate 
that  the  rate  of  irreversible  entropy  production  is  a result  of  the  various  dis- 
sipative phenomena:  the  first  term  is  the  contribution  due  to  mass  transfer. 


Suppose  that  the  chemical  kinetics  of  the  reacting  mixtures  may  be 
described  in  terms  of  a set  of  a.,.  chemical  reactlonc,  which  may  be 
written  symbolically  thus: 


A1S[1]  + Aft  [2]  +•*•  — e>Jl]  + bZi[2l  +• 


(Uta) 


in  which 


indicates  the  chemical  formula  for  the  I1*1  species  and  A 


th  * i* 

and  B}8  are  the  iutegrai  Btoicluometric  factors  in  the  »w  reaction  fox  the 
i1*1  substance.  Then  if  kB  and  kj,  are  the  rate  constants  for  the 


forward  and  backward  reactions  respectively, 
r i.  Lht*  \ ( \.  ' £ r 


TJid  Y8  are  defined  by: 


(M,  C-)  - (kit, 6,1  ft 


Y,  « - Z (6,,-Ato)^, 

!»■! 


(llZt) 

(1-1  Zc) 


in  which  the  fj  are  fugacities. 

The  quantities  r8  and  Y8  are  then  •’fluxes'1  and  "driving  forces" 
respectively. 
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t he  £4Cond  due  to  viscous  effects,  the  third  die  to  heat  transfer,  and  the  last 
due  to  chemical  reaction.  It  should  be  noticed  that  according  to  E<q.  (1.  12) 
the  rate  of  entropy  production  g can  be  written  as  the  sum  of  products 
of  fluxes  and  driving  forces  (the  latter  are  lametimeB  called  ♦'affinities*' ) . 

This  is  an  important  result  which  can  be  used  in  conjunction  with  the  thermo- 
dynamics of  irreversible  processes  to  obtain  information  concerning  the 
relation  between  the  fluxes  and  the  transport  coefficients.  This  is  discussed 
further  in  §2. 

C A ppUc  ability  of  the  Equations  0/  Chan  fig 

Although  the  equations  of  change  as  given  above  are  general  and  apply  to 
any  fluid,  they  .~ra  clearly  useful  only  under  Buch  conditions  that  it  is  physical- 
ly meaningful  to  speak  about  point  properties.  The  definitions  of  local  density, 
velocity,  and  temperature  arc  reasonable  only  if  the  fluid  can  be  regarded  as 
behaving  as  a continuum.  When  there  are  large  differences  in  the  macroscopic 
quantities  over  distances  of  the  older  of  a mean  free  path,  the  distribution  of 
the  velocities  of  the  molecules  deviates  considerably  from  a Maxwellian 
distribution.  Under  these  conditions,  the  concepts  of  local  composition,  velocity, 
and  temperature  become  meaningless.  Two  well-known  examples  of  a fluid  under 
Btii’n  a condition  are:  (i)  an  extremely  dilute  gas  (or  "Knudsen  gas")  in  which  the 
dime  as  id  ns  of  the  containing  vessel  or  an  object  immersed  in  the  gas  Are  of  the 
same  order  of  magnitude  as  the  mean  free  path,  and  (ii)  a shock  wave  in  wh  :h 
the  macroscopic  properties  as  a function  of  the  distance  undergo  an  abrupt  change 
within  a distance  of  a few  mean  flee  paths.  In  both  these  c*Ben  it  is  meaningless 
to  speak  of  point  values  of  the  macroscopic  variables. 
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One  might  well  question  the  applicability  of  the  equations  of  change  to  the 
description  of  turbulent  flow.  As  the  flow  conditions  change  ffOMn  laminar  ts 
turbulent,  small  eddies  appear,  and  there  is  superimposed  on  the  overall 
flow  a complex  pattern  of  more  or  less  random  motions.  The  dimension* 
of  these  eddies  are  always  large  compared  to  the  mean  tree  pftth,  *o  that 
the  turbulent  motion  is  macroscopic  rather  than  molecular.  Accordingly, 
the  concept  of  the  fluid  as  a continuum  remains  valid,  and  the  ordinary 
equations  of  change  may  be  used.  In  these  equations  the  variables  refer  to 
instantaneous  values  at  a point.  For  most  practical  purposes,  however,  we 
are  interested  in  the  values  of  these  quantities  averaged  pver  a time  long 
compared  with  the  period  of  fluctuations.  Hence,  in  turbulent  flow,  the 
equations  of  change  given  above  sire  modified  by  time -aver aging  to  obtain 
relations  among  the  average  macroscopic  variable*.  When  thi*  is  done 
one  obtains  a set  of  time - averaged  equations  of  change  which  dilfer  in.  two 
respects  from  the  original  equations:  (1)  the  variables  and  fluxes  appearing 
ir.  the  equations  are  now  time -averaged  quantities-  »ad  {ii)  there  arc  addi- 
tional terms  r.nsociatcd  with  the  correlation*  of  the  tluctuation*  in  the  various 


physical  properties  (in  the  equation  of  motion,  for  example,  one  obtain*  the 
"Reynolds  stresses").  The  explaining  of  tlwse  added  terms  ia  one  of  the 
principal  problems  of  turbulence.  The  simplest  theories  for  this  are  the 
mixing  length  theories,  which  are  analogous  to  the  mean  free  path  theories 
of  elementary  kindle  theories.  Another  even  simpler  approach  is  to  presume 
vhai  turbulent  transport  is  analogous  to  laminar  transport  and  thus  to  introduce 
(.he  gad/  coefficients  cf  diffusion,  viscosity,  and  therms)  conductivity.  These 
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quaatltles  are  not  constant*  of  the  fluid,  but  depend  on  the  position  in  the  fluid, 
the  Reynr.lds  number  for  the  flow,  and  also  on  the  scale  of  the  turbulence.  Hence 
the  use  of  eddy  transport  coefficients  is  not  too  satisfactory.  A summary  of 
the  present  state  of  experimental  and  theoretical  knowledge  of  the  Reynolds 
stresses  and  eddy  viscosity  has  been  prepared  by  Schlichting  (35);  and  the  sub- 
ject of  eddy  diffusion  is  summarised  by  Sherwood  and  Pigford  (39).  Other 
recent  work  on  mass  and  heat  transfer  in  turbulent  systems  has  been  done  at 
California  Institute  of  Technology  by  Schllnger,  Sage,  and  collaborators  (8, 

36.  37). 

In  flow  systems  consisting  cf  more  than  one  phase,  such  as  one  encounters 
in  industrial  contact!*  g processes,  the  equations  of  change  are  valid  in  each  of 
the  phases  concerned,  but  the  solutions  to  the  equations  of  change  have  to  Mmatch 
up"  at  the  interfaces  between  the  phases.  At  the  present,  there  is  only  merger 
knowledge  as  to  ho**  the  solutions  to  the  t ^uations  of  change  should  be  joined 
at  the  phase  boundary^  inasmuch  as  very  little  is  known  about  the  mechanism 
of  in.terpha.se  mass,  momentum,  and  energy  transfer. 

Recently  two  groups  o t investigators  have  studied  the  problem  of  the 
mechanism  of  interphase  mass  transfer.  Scott,  Tung,  and  Drickamer  (38) 
solved  the  differential  equations  for  equimolal  diffusion  across  in  interface, 
taking  into  Account  the  fact  that  the  interface  itself  would  provide  an  additional 
resistance  to  the  diffusions!  process . Subsequent  radio-tracer  experiments 
of  Tung  and  Drickamer  (40)  indicated  that  the  interfacial  resistance  can  in- 
fluence very  markedly  the  concentration  profiles  during  diffusion.  Another 
approach  to  the  problem  was  taken  by  Emmert  and  Pigford,  (13)  who  studied 
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absorptian  in  wetted  wall  towers  in  terms  of  accommodation  coefficients. 
They  also  concluded  that  interfacial  resistance  is  of  considerable  importance. 
Certainly  this  question  of  the  mechanism  of  interphase  mass -transfer  needs 
further  study.  One  aspect  of  the  subject  of  interphase  mass-transfer  has 
been  exhaustively  studied  by  chemical  engineers,  namely  that  of  correlating 
overall  mass -transfe  r data  with  flow  variables  according  to  dimensional 
considerations.  An  excellent  summary  and  evaluation  of  this  material  has 
been  given  b*  Sherwood  and  Pigford.  (39) 
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§2  The  Fluxes  in  Terms  of  the  Transport  Coefficients 

The  equations  of  change  in  terms  of  the  flux,  vectors  as  given  in  §)  have 
been  derived  by  means  of  non.  equilibrium  statistical  mechanics.  According 
to  statistical  mechanics  the  flux  vectors  in  turn  can  be  expressed  in  terms  of 
integrals  which  involve  distribution  functions.  Ultimately  one  would  like  to 
follow  through  this  line  of  attack  and  obtain  expreasu  ns  for  the  distribution 
functions  and  thereby  express  the  flux  vectors  in  terms  of  the  forces  between 
the  molecules  In  the  fluid  and  the  various  derivatives  of  the  macroscopic 
properties.  The  complete  development  of  the  connection  between  the  flux 
vectors  and  the  inter  mole  rular  forces  has  been  worked  cut  only  for  certain 
special  types  of  molecules  and  for  certain  limited  conditions;  (i)  For  the  dilute 
monatomic  gas  the  Chapman-llnskog  theory  (6)  gives  the  flux  vectors  explicitly 
in  terms  of  the  transport  coefficients,  which  in  turn  are  expressed  in  terms 
of  intermolecular  forces  as  described  in  |j4.  A considerable  amount  of 
computational  work  has  been  done  for  the  dilute  milticomponent  monatomic 
gas  mixtures  with  small  gradients,  and  it  is  possible  to  predict  quite  accurate 
ly  the  flow  properties  under  both  classical  and  quantum  conditions,  (ii)  For 
the  dilute  polyatomic  gas  a general  theory  has  been  developed  by  Wang  Chang 
and  Uhlenbeck  (42)  and  also  by  de  Boer  (10)  but  no  numerical  results  have  yet 
been  obtained  whereby  one  can  relate  the  transport  coefficients  to  the  infcor- 
molccular  forces.  ( iii)  For  the  dense  monatomic  gas  formal  results  have 
been  obtained  by  Kirkwood  and  coworkers  (21)  ai:d  by  Born  and  Green  (3) 
but  much  work  needs  to  be  done  in  order  to  make  these  results  of  practical 


value . 
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For  a fluid  made  up  of  any  kind  of  molecules  an*l  under  any  arbitrary 
conditions  of  temperature  and  pressure,  it  is  not  possible  at  the  present  time 
to  use  statistical  mechanics  to  obtain  the  form  of  the  flux  vectors  and  to  derive 
expressions  for  the  transport  properties  in  i;rms  of  inter  molecular  forces. 

For  suen  a “general  -fluid",  however,  a certain  amount  of  information  may  be 
obtained  from  the  thermodynamics  of  irreversible  processes.  It  is  possible 
to  derive  the  form  of  the  flux  vectors,  and  by  means  of  the  OnBager  "reciprocal 
relations"  to  obtain  certain  relationships  between  the  transport  coefficients. 

It  should  be  emphasized,  however,  that  the  thermodynamic  approach  cannot 
lead  to  any  expressions  for  the  transport  coefficients  in  terms  of  the  molecular 
properties.  Hence  we  summarize  In  this  section  the  most  important  useful 
results  which  are  obtained  about  the  flux  vectors  from  the  thermodynamics 
of  irreversible  processes  and  other  sources. 

a.  The  Thei modynamics  of  Irreversible  Processes 

In  Eq.  (1.  11)  it  was  found  that  the  rate  of  entropy  production  could  be 

written  as  a sum  of  products  of  fluxes  and  driving  forces.  That  is: 

3nt«-  *1|  » 12 

Cl  = ( 1 / T ) 2Z,  Jn-X-r,  (2.1) 

rul 

in  which  the  Jn  are  the  various  fluxes  and  the  Xn  are  the  various  driving 
forceB  or  "affinities.*'  These  fluxes  and  forces  are  summarized  in  Table  I. 

In  general,  any  driving  force  can  give  rise  to  any  of  the  fluxeB.  For  exampl 
if  th^re  is  in  a syuteei  both  & gradient  and  a temperature  gradient,  the  fol- 


lowing fluxes  are  observed  (i)  mass  fliix  due  to  the  concentration  gradient 
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Table  I 
Fluxes  and  Forces  in  < 

i Flawing  Fluid 

Flux.  Jn 

Force,  X 

n 

[ ‘ 

it 

i 

l> 

r*  - 

r 

- Vi* 

-VI  n T 

rs 

*Ys 

Notes:  The  quantity  l\  . is  defined  in  Eq.  (1.13),  and  r8  and  Y# 

are  defined  in  Eqs.  (1.  12b).  The  subscript  i goes  from  1 
to  nc  (the  total  numb*.’  of  chemical  species  present),  and  a 
goes  from  I to  nr  (the  total  number  of  chemical  reactions 
taking  place).  The  quantities  r8  and  YB  are  scalars; 
and  and  »»»nriat»H  fnrr.es  are  vector  quantities 

t . 

with  three  components  each;  and  p'v'  and  V*  are  second 

a 

order  tensors  with  nine  components  each.  Hence,  altogether 
we  arc  concerned  with  3n(.  + ny  + 12  flukes  and  an  aqua.1 
number  of  forces. 


WIS-ONR-7 
12  Miy  1954 
-15- 


(ordinarv  diffusion);  (ii)  energy  flux  due  to  the  temperature  gradient  (thermal 
conductivity);  (ii4)  mass  flux  due  to  the  temperature  gradient  (thermal  dif- 
fusion the  "Soret  effect");  and  (iv)  energy  flux  due  to  the  mass  gr  «:  client 
(the  diifusion-ther moeffect  o.’  the  "Dufour  effect").  We  see  thus  that  there 
are  two  types  c#f  effects- --direct  effects,  such  aa  (i)  and  (ii),  and  coupled 
effects,  such  as  (iii)  and  (iv). 

It  is  assumed  as  a postulate  of  the  thermodynamics  of  irreversible 
processes  that,  for  situations  not  too  far  removed  from  equilibrium,  the 
fluxes  may  be  written  in  the  form: 


(ZZ) 


which  states  that  the  fluxes  axe  linear  functions  of  the  driving  forces.  The 

anm  are  referrr  d to  as  the  phenomenological,  coefficients.  The  diagonal 

elements  apn  are  the  coefficients  which  represent  the  direct  effects;  the 

off -diagonal  elemenets  anm  (n  ± m)  are  the'phenomenological  cocfficientsJ 

for  the  coupled  effects.  The  expressions  for  the  fluxes  obtained  by  statistical 

mechanical  considerations  are  of  this  form.  Also,  such  "linear"  relations 

are  found  in  connection  with  oiher  types  of  physical  and  chemical  phenomena— 

for  example,  in  thermo  magnetic  effects,  thermo  galvanic  effects,  electrokinetic 
II 

effects,  etc.1  Linear  relations  should  always  be  adequate  when  the  systems 
are  sufficiently  close  to  equilibrium.  Hence  the  postulate  of  Eq.  (2.  2)  defines 
the  range  of  applicability  of  the  thermodynamics  of  irreversible  processes. 


Complete  discussions  of  these  applications  inav  be  found  in  the  recent  mono- 
graph of  de  Grant  (12). 


WIS-ONR-7 
1?.  May  1954 
-16- 


The  fundamental  theorem  of  the  thermodynamics  of  irreversible  processes 
Is  due  to  Onsager  (31).  This  theorem  states  that  for  & ''proper  choice"  of  the 
fluxes  and  forces 


a 


wn 


l'? 


These  equations  are  comnaorly  referred  to  as  'Onsager*s  reciprocal  relatione 

and  are  useful  for  relating  coupled  effects  such  as  the  Soret  and  Dufonr  effects. 

Let  us  now  consider  the  application  of  the  linear  law  of  Eq.  (2.  2)  and  the 

Onsager  reciprocal  relations  to  the  study  of  transport  phenomena  in  fluids, 

for  which  a proper  choice  of  the  fluxes  and  forceB  is  indicated  in  Table  1. 

When  the  linear  Jaw  for  the  fluxes  is  applied  to  the  flux  vectors  which  are 

equivalent  to  the  first  approximation  in  the  Chapman-Enskog  kinetic  theory. 

That  is,  one  obtains  expressions  for  the  flux  vectors  which  are  linear  in  the 

first  derivatives  of  the  macroscopic  variables  and  do  not  contain  higher 

derivatives  or  powers  of  the  first  derivatives. 

Hence,  if  the  system  under  consideration  in  in  a state  not  too  far  removed 

from  equilibrium,  the  linear  law  of  Eq.  (2.  2)  indicates  that  each  of  the  components 

of  the  fluxes  in  Table  I may  be  written  ae  a linear  combination  of  all  ti.e  components 

of  all  the  forces  in  Table  I.  However,  if  the  system  is  isotropic,  it  raay  be  shown 

that  those  terms  which  correspond  to  a coupling  of  tensors  whose  orders  differ 

by  an  odd  number  do  not  occur.  Thue  no  coupling  occurs  between  jj  or 

(both  vectors  or  first-order  tensors)  and  p^v)  (a  second  order  tensor)  or  r0 

(a  scalar  or  zero-order  tensor).  Coupling  does  occur  between  the  heat  flux 
ih) 

er  ' and  the  mass  flux  . resulting  in  the  Soret  effect  and  the  Dofour  effect. 
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Coupling  Also  cccuxe  between  the  momentum  flux  anti  the  chemical  reaction 
rates,  since  there  is  a difference  ol  two  in  the  order  of  the  tensors.  No 
experimental  observations  of  such  an  effect  have,  however,  been  reported, 
b.  Momentum  Flux  and  the  Viscosrtv  Coefficients 


The  application  of  the  liiu-«...r  law  (Eq.  (2.  2))  to  the  components  of  the 

viscous  pressure  tensor  p(v)  and  the  force  - Vv  gives; 

2 * 

V 1 UL 

rhk  = - A A,  (do>-k/6i<K)  (2.4) 

n«i  m-i 

hk 

The  set  of  phenomenological  coefficients  anm  is  a matrix  of  eighty-one 
elements,  but  all  of  Ihpee  are  not  independent.  From  Qnso^er's  relations 
one  sees  that  aj^  * 41k1  ; and  since  the  pressure  tensor  is  defined,  to  be 

symmetric,  it  can  be  el- iwn  that  = »nm  . Then  Eq.  (2.  4)  can  bo 

simplified  to: 


bf  - -mlj,  [(ivK/bxh)  + (H/d>-k)] 


(2.5) 


n»i  m«i 


sp  that  there  are  only  thirty-s^x  independent  ‘COfff  \<3  jjen+*  ^ 

If  the  further  assumption  is  made  that  the  fluid  is  inotropic,  then  It  cjin 

hk 

shown  that  the  number  of  independent  anm  is  reduced  to  two.  The  pressure 
tensor  then  has  the  form: 

j>IV’  - -P-  [(Vs)  + (Vur)1  ] + ( - K)(V.jr)  r 


lik 

anm  • 


in  which  |A-  and  K»  are  the  two  independent  phenomenological  coefficients, 
which  are  called  the  coefficients  of  shear  and  bulk  viscosity  respectively. 


f 


Typical  diagonal  and  non-diagonal  elements  of  this  tensor  are; 


h " 

« XX. 


'i  x~  l'''! 

.3*  A bx 


$ , 1B1 


?;£  / .< 


k'V> 


* -J&  + 

rv  oy 


*5i) 
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The  coefficient  of  shear  viscosity  is  important  in  fkiwo  is  which  successive 
plane  layer  a of  the  fluid  move  with  different  velocities.  The  coefficient  of 
bulk  viscosity  is  important  in  the  pure  expansion  of  A fluid. 

In  a dilate  a fc  made  up  of  molecules  without  f *•  * — r rsl  wtfjjSwwStnf  fsrcicdmn, 
the  coefficient  of  bulk  viscosity  is  identically  Kero.  It  *h*wn  that  the 

presence  of  internal  degrees  of  freedom,  even  in  a dilute  gas,  iitttVduca*  a 
finite  bulk',  viscosity  t and  the  bulk  viscosity  has  also  !*•*«.  to  bs  nan-r-erc 

in  a dense  gas  or  liquin.  Nevertheless  the  effect!*  AdtfeXiy  mwt&Xx,  fluid  it  is 
ofte&nfiuffic iently  accurate  to  neglect  the  bulk  viscosity.  A survey  of  ttur  topic 
of  bulk  viscosity  has  been  given  by  Karim  and  Ro  8 Sinks  Adi  {S4)i  tf  ihs  fluid 
under  consideration  is  incompressible  (_then  ( V«y)  »0),  tka  last  W m of 
Eq.  (Z.6)  vanishes,  and  the  coefficient  of  bulk  vissoaity~e»ooiB  uOfrotttog  into 
ths  equation. 

For  the  flow  of  a fluid  in  the  x-dircction  with  a velocity  gradient  in  the 
y -direction  only,  Eq.  (Z.o)  simplifies  to: 


= -p.  <bv,lby)  (Zl) 

W | 

This  is  the  well-known  Newton’  8 law  of  viscosity  which  *tAt#A  that  the 
resistance  to  viscous  force3  is  proportional  to  the  vcU>Uy  grAdtettk  the 
constant  of  proportionality  being  the  coefficient  of  shsijjp  *4»sssity.  The  co- 
efficient of  shear  viscosity  is  determined  experimentally  by  considering  the 
flow  of  a fluid  in  a system  of  aulficlently  simple  geometry  that  the  equation 
of  motion  can  be  solved  after  the  expression  for  ihe  vice  cos  sirens  given  in 
Eq.  (2  7)  lias  b**en  substituted  into  if.  In  Ibis  way  one  analyses  thq  data  for 
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capillary  tube  viscosimeters  in  terms  of  Poisemll*'*  law,  and  the  data  for 
concentric  cylinder  (Couette)  viscosimeters  in  terms  of  the  Margules 
equation. Eq.  (2.7)  then  forms  the  basis  for  the  entire  subject  of  viscosimetry 
of  newtonian  fluids. 

Thdcfc  axv  many  fluid;  which  do  not  behave  according  to  the  above 
f'newtonian*4)  form  of  the  pressure  tensor,  and  such  fluids  are  referred  to  as 
"non-newtonian".  Examples  of  non-newtotdan  substances  are  glass,  plastics 
melts.,  drilling  muds,  '’bouncing  putty,*4  toothpaste,  and  concrete.  For  these 
substances  the  assumption  of  small  gradients,  made  in  connection  with  Eq. 

(2.2),  is  net  tenable.  Hence  tlie  pressure  tensor  depends  upon  higher  deriva- 
tives of  the  velocity  components  and  on  higher  powers  of  the  velocity  gradients. 
The  inclusion  of  these  additional  terms  would  require  the  introduction  of  addi- 
tional phenomenological  coefficients;  approaches  of  thin  sort  have  been  quite 
limited. 

In  connection  with  the  establishment  of  this  relation  between  the  shear 
stress  and  the  shear  rate  for  non-newtonian  flow,  considerable  success  has 
been  obtained  bv  Eyring  and  his  collaborators  (15)  who  have  applied  the  meth- 
ods of  reaction  rate  theory  to  various  postulated  physical  pictures  of  the  elemen- 
tary processes  going  on  Inside  a fluid.  It  is  this  theory  which  predicts  that  the 
shear  viscosity  depends  on  the  external  force  according  the  the  hyperbolic  sine 
law. 

Actually  moot  of  the  studies  of  the  flow  behavior  of  non-nevtcnian  substances 
have  been  baaed  on  simple  analytical  assumptions  concerning  the  relationship 
between  the  sheai  stress  (p^)  and  the  shear  rale  (dvx/5y).  This  ie  also 
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tantamount  to  Baying  that  come  assumption  is  made  concerning  the  depenaence 
o£  the  coefficient  of  viscosity  on  the  impressed  shear  force.  By  means  trf  these 
assumed  "models",  one  can  solve  the  equation  of  motion  for  various  simple  ge- 
ometries (such  as  capillary  tubes,  concentric  cylinders,  etc.)  and  obtain  for- 
mulae which  are  useful  for  the  viscosimetry  of  non-newtonian  substances.  (33). 
Formulae  of  this  sort  are  available  for  various  idealized  non-newtonisn  medels, 
such  as  the  Bingham  plastic,  thixotropic  substances,  and  power-law  models 
(which  just  assume  that,  the  shear  rate  varies  as  some  power  of  the  shear  stress, 
nr  vice  versa).  Unfortunately  many  investigators  continue  to  analyze  the  flow 
data  for  non-newtonian  substances  in  terms  of  PoiseuMle's  law  and  other  new- 
tonian  formulae,  thereby  obtaining  "effective  viscosities'*  which  are  valuable 
only  for  analyzing  flow  in  systems  geometrically  similar  to  the  viscosity  meas- 
urement apparatus.  The  study  of  the  flow  and  deformation  of  these  non-newtonian 
substances  is  now  generally  classified  under  the  heading  of  "rheology* " This 
subject  includes  the  behavior  of  substances  under  stress,  all  the  way  from  pure 
elastic  deformation  (Hooke's  law  to  pure  viscous  flow  (Newton'S  law). 

c . Mass  and  Heat  Flux  and  the  Coefficients  of  Diffusion  and  Thermal 
C onductivlty 

The  application  of  the  linear  law  (Eq.(2.  2))  to  the  components  of  the  heat 
and  mass  flux  vectors  and  the  corresponding  driving  forces  given  in  Table  I 
gives  (after  some  algebraic  rearrangement): 


nt 


_(b)  ..  r-y  1 t 

3 = - V In  T 

- P - a0j  ’?j)  4j 

J_i  =•  — 3(o  V In  T 

- F>  2 (atj  /j>j>  dj 

(>•£> 
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in  which  the  driving  forces  d.  , closely  related  to  the  .A  ^ , are  defined  thu#,; 

dj  - (iV'p)  Aj  - (fj/p)  V hp  + (f-j/fp)  Zpkfk 


>v 


+ UfjVj  /mj)  - (pj  /p)  ] V In  p 

- Cpj/ppXffj  - Zpkfk) 


k„  r-C  ' (CIO) 

In  the  above  equations  the  index  0 refers  to  the  temperature  variable;  the  other 
the  other  subscripts  (i,  j,  k)  refer  to  the  chemical  specie^  present  in  the  mix- 
ture. The  following  restrictions  have  to  be  placed  gn  the  a^j:  (f)  because  of 
the  Onsager  relations  the  ay  are  symmetric;  (ii)  because  of  tb*  fact  that  the 
j j are  not  all  independent  (i.e.,  *?’.]  = G ),  the  ay  must  in  addition  satis- 

V'  XT' 

fy  the  relations:  Z.{ ay  r j - 0. 

The  app-ic».«ion  of  the  Onsager  relations  to  the  phenomenological  coefficients 
wlf.J  one  0 subscript  gives; 


ai0  - <3oi  31  Di  (?-Vi) 

and  -we  have  equated  these  elements  to  the  'multicomponent  thermal  diffusion 
coefficient,  ••  D*^  . Thus  the  Reciprocal  relation  in  Eq.  (2.  11)  gives  the  con- 
nection between  the  coupled  effects  known  as  thermal -diffusion  and  the  diffusion- 
thermo  effects. 

The  application  of  the  Onsager  relations  to  the  remaining  phenomenological 
coefficients  gives: 


U ¥=  0 

J 
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We  now  have  to  relate  theBe  quantities  to  the  "smjlti««>**pe9u»Qi  diffusion  coef- 
ficient*. M D^j.  but  this  has  to  he  done  in  auch  a -vfay .the*  Eq*  {Zi.  iZ't  ia  obeyed. 
Inasmuch  as  the  Dj  j are  customarily  defined  so  that  the  Q|  ^ 4^  identically 

aero  (9).  the.  relationships  between  the  and  the  a^j  are  ir Other  involved: 

rv 


aij  = (n^/p'b)  Djj  + Z 


nc 


(213) 

(2.14) 


Dii  = -(fp/hl^rY,j)[(l/ft)atj  + (</ft)Zaik 

J * k*>1 

Mi 

It  should  be  noted  that  these  Dy  are  not  symmetric  with  respect  to  aa  inter- 
change of  the  subscripts  i and  j . 


For  ateo-compcnent  mixture  the  “mvUticomponeot  diffusion  coefficients” 
Bjj  become  the  ordinary  binary  diffusion  coefficients  Sbij  , and  then 
c&ij  =<£)ji  • For  two  -component  systems  it  is  alee  customary  to  define  a 
Mth*rmal  diffusion  ratio,  !?  kj,  according  to  the  relationt 

,i  'it  |M  ’9I  i 

Ic_  » (j)/n2m,mt)(  D(T/£)I2.)  c - (f  /*ilrn;rrn)(I>t  /&it)  (2.14a) 


For  this  definition  of  k^,,  component  goes  to  the  cold  re  clot  If  k.j.:>  0 
and  to  the  hot  region  if  k-p < 0.  The  diffusion  coefficients  obi j and  the 
thermal  diffusion  ratio  k^,  defined  thus  are  in  agree  gaumt  with  the  ac- 
cepted definitions  of  the  previous  -worker**,  and  in  particular  with  the  text- 
books of  Chapman  and  Cowling^aad  Grew  and  IbbsfiQ. 

For  a three -component  mixture  the  multicomponent  diffusion  coefficients 
are  not  equal  to  the  ordinary  binary  diffusion  coefficients.  For  example, 
the  kinetic  theory  of  dilute  multicomponent  gas  rstixture^'kho w ft  that; 


Diz.  ~ 


1 1 + 

L 


ru  [(  r'r'o/m7.)  §)15  — 3 1 


(2.14b) 
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We  now  can  substitute  the  results  of  the  application  of  the  Onsager  relations 
(Eq.  (2.  11)  and  (2.  13))  into  the  expressions  for  the  near  and  mass  fluxes  (Eq. 
(2.  8)  and  (2.9))  to  obtain: 

f1-  -aooVloT  - f>Z(D,T/pf)dj  (z-<5> 

U = -Di'VInT  t (nl/f)  f,  m^j  D(j  dj  (2.16) 

in  which  the  phenomenological  coefficient*  aQ  q still  appears.  The  quontity  <^0o/'l' 
is  not  quite  the  coefficient  of  thermal  conductivity  as  customarily  dofined,  as 
can  be  seen  in  the  following  way:  Imagine  that  Eq.  (2.  16)  is  solved  for  the  d t 

and  that  these  quantities  are  in  turn  substituted  into  Eq.  (2.  15).  This  will  then 
indicate  that  has  other  terms  proportional  to  VT  besides  the  term  con- 

taining a . When  all  of  these  terms  a.-e  combined,  the  coefficient  of  the  re- 
0 0 

suiting  term  is  detined  to  be  the  coefficient  of  thermal  conductivity.  This  par- 
ticular choice  for  the  definition  of  is  made  for  the  sake  of  consistency  with 
existing  treatises  on  kinetic  theory,  such  as  that  of  Chapman  and  Cowling  (6)  . 

When  the  rearrangements  just  described  have  been  performed,  one  finally 
obtains  for  the  components  of  the  energy  flux  vector  (as  given  in  Eq.  (1.6)) 
the  following: 


q,n  = -'A  VT  (2.17) 

r< 


I*1 

(T) 

The  first  term,  q , ie  the  contribution  to  the  energy  flux  due  to  thermal 


conduction  and  is  the  coefficient  of  thermal  conductivity.  The  tu-cond  term, 
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(d) 

q , l*  the  energy  flux  associated  with  the  diffusional  processes  occurring. 
Thacs  tws  toroM  ur<e  the  most  important  contributions  to  the  energy  flux. 

The  third  to? ^ ii  related  to  {but  not  equal  to). the  energy  flux  resulting 

from  the  JDufo^r  effect.  Actually  Eq.  (2.19)  is  rigorously  applicable  only  to 
a dilute  gas.  mixture.  .The  expression  for  far  a general  fluid  is 

extremely  complex  and  the  effect  is  small.  ui 

In  the  above  set  of  equations  ve  have  given  no  expression  for  , the 

energy  flust  d*#  to  Mdintion. This  set  of  aquations  applies  only  to  systems 
near  equilibrium*  1£  radiative  energy  transport  is  taking  place  and  if  the 
radiation  is.  in  equilibrium  with  the  matter,  then  the  expression  Lor  <^r) 
is  of  the  same  form  as  Eq.  (2.  17)  for  q^'  . This  condition  ic  usually  not 
a slip  fled  and  radiative  energy  flux  depends  in  a detailed  manner  on  the 

£ 

frequency  amd  i^teoaity  of  the  radiation  and  the  properties  of  the  material.'3 

From  Eq*,.  (2^16)  and  (2,  10)  one  obtains  the  virion*  contributions  to 
the  mass  flux  vector(«a  given  in  Eq.  (1.4)),  thus; 

«jt  - t (*Vjp)  "Wj D<j  [(n}''p) Z (tyj  Vp*  J (2.2o) 


fa  _ >> 

J>  ” +SW  Z "W*>j  I^j  [(^Vj  /m,)  - fPj  /f )]  V In  f> 


/e  • A 


hi  he 

]i{r)  « * (HVf)  Z D.j  [fp:  /f>p)(p  fj  ” Z ) J (2^*22) 


i r 

iv 


(T) 


Di  V in  T 


(Z.25) 


.1  The  subject  of  radiative  heat  transfer  is  discussed  further  in  MTGL,  § l i . 3 . 
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The  first  contribution,  , is  the  most  important  one--the  flux  due  to  a 

i 

concentration  gradient.  This  expression  indicates  that  the  correct  driving 

force  for  diffusion  is  a rather  complicated  function  of  the  chemical  potentials 

of  all  the  components  in  the  mixture.  The  mass  flux  due  to  a gradient  in  the 

static  pressure,  . has  been  studied  only  slightly  because  of  ftp  rela- 

tively  negligible  importance.  The  mass-flux  due  to  forced  diffusion,  , 

is  important  when  systems  containing  ions  are  under  consideration.  And  the 
(T) 

last  term,  j:  , is  the  contribution  to  the  mass  flux  due  to  thermal  diffusion. 

This  contribution  can  be  neglected  except  under  conditions  of  large  terpera 
ture  gradients  such  as  are  carefully  maintained  in  devices  such  as  the  Cluuiu 
Dickel  'column  (16). 

For  two-component  systems  the  mass -flux  vectors  become  considerably 
simpler.  The  general  expression  to~r  the  flux  of  component  "1M  iSr 

ji  - » - K/j>)  (ds  + kT  V lr,  T)  fe.fc*.') 

in  which  we  have  also  indicated  the  relation  between  the  flux  vector  and  the 
diffusion  velocity,  . The  latter  is  the  velocity  with  which  component 

" 1"  moves  with  respect  to  the  mass  average  velocity.  Inasmuch  -e  the  bi 

\ ' 

nary  diffusion  constants  are  symmetric  ( c an^  ®ince  ^ 


it  follows  that; 
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Y2  -Y,(  * (nVn!  n2.)  §J|?  (4,4-  ^ l,r',^Y  (2.25)f 

This  equation  is  the  usual  starting  point  for  discussion  of  thermal  diffusion. 

For  a two-  component  system  in  which  thermal,  pressure,  and  forced  dif 
fusion  effects  can  be  neglected,  the  expression  for  the  mass  flux  given  in 
Eq.(2.24)  becomes; 

j,  = ~(rrn£/j5f)m,m2^iZ  V|2,  = -(nVp)(n£kT/p)rr>,rn^5bu  (ZU) 

in  which  aj  is  the  activity  of  component  ”1".  From  this  it  can  be  seen  that 
it  is  the  chemical  potential  or  the  activity  which  should  correctly  be  considered 
to  be  the  driving  force  for  two-component  diffusion.  When  the  additional  asnmp 
tion  is  made  that  the  components  are  dilute  gases,  then  Eq.(2.  26)  becomes; 

ji  ■ - (nVp)  V*i  (2.2.7) 

This  expression  can  be  used  in  systems  where  there  are  temperature  gradients, 
provided  that  they  are  not  sufficiently  large  that  the  thermal  diffusion  contri- 
bution needs  to  be  added. 


For  a dilute  monatomic  gas  this  result  has  been  generalised  for  multicom- 
ponent mixtures: 

2 Ov*)  • 3>*J  X Y/A)  -YtA>)  - (Z'25^ 

3-'  ru 


di  ~ (v  In  r)  2 (n.nj/ n^ij)  [(Dj/Vijrnj)  - (Dt/njmij] 

3 

For  systems  at  constant  temperature  and  pressure  and  with  no  external 

forces,  this  simplifies  to; 

IV 


- Vpt  - 2 oi,j  ci  Cj  (Kw-  If) 


(2., 25  b) 


1 W| 

in  which  the  "resistance  factors,  " of  , are  defined  as  (RT)Z  /p56;  .,  anil 
Ci  is  the  molar  concentration  of  epeM.es  i . This  result,  originally  due 
to  Maxwell  and  Stefan,  is  commonly  used  as  the  stariiug  point;  far  discussions 
of  multicomponent  diffusion  (203  • 
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It  has  been  previously  mentioned  that  the  flux  vectors  defined  abeve  are 
fluxes  with  respect  to  a coordinate  system  moving  with  ths  mass  average 
velocity  of  the  gas.  Sometimes  it  if  convenient  to  define  fluxes  with  respect 
to  other  coordinate  frames.  Two  examples  will  illustrate  the  importance  of 
specifying  the  reference  frame  for  the  fluxes;  (i)  The  flux,  in  moles 

per  unit  area  n*i  unit  time  vith  respect  to  the  molar  average  velocity  (under 
the  same  restrictions  of  Eq.  (2.  27))  is; 

J,  = - (paWRT)Vx,  (ZZH) 

(ii)  The  molar  flux  of  component  " 1M  with  respect  to  a stationary  coordi- 
nate system  is; 

N,  m — (p&iz/Rf)  Vx|  +•  (bl|  +*  (2.29) 

The  first  term  on  the  right  hand  aide  is  the  molar  flux  of  component  •>!"  by 
diffusion,  and  the  second  term  is  that  due  to  the  muss  motion  of  the  fluid. 
Eq.(2.29)  is  the  basic  starting  point  for  the  diBcusBion  of  chemical  engineer- 
ing mass-transfer  problems.  A large  number  of  the  Important  applications 
involve  either  of  two  assumptions;  (i)  diffusion  of  *' 1”  through  stagnant 
M2M  , so  that  Nz  = * '•  anti  (**)  «quimolal  coumerdiffusicnfor  which  = 

-N^  - The  applications  of  Eq.(2.29)  to  various  important  industrial  diffu- 
sional  operations  have  been  given  by  Sherwood  and  Pigford.(3^)). 

The  expressions  given  in  this  section  for  the  mo  men. turn  flux  (Eq.(2.  6))^ 
tne  components  of  the  heat  flux  (E<>s.(2.  1Y)  to  (?.  19)),  »n«’L  the  components 
of  the  mass  flux  (Eqs.(2.  20)  to  (2.  23)  ) are  general  results  obtained  from  the 
thermodynamics  of  irreversible  processes  and  which  apply  to  both  gases  and 
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liquids.  They  are,  of  cours;  consistent  with  the  expressions  for  the  flux 
vectors  for  dilute  monatomic:  gases,  which  have  been  previously  obtained 
from  the  rigorous  kinetic  theory  of  gases. 

When  the  expressions  for  the  components  of  the  mass  flux  (Eqs.(2.  20) 

to  (2.  23)  ) are  substituted  into  the  equation  of  continuity  for  component  i 

(Eq.(l.  la)),  one  obtains  the  diffusion  equation  in  its  most  general  form. 

When  it  is  assumed  that  thermal,  pressure,  and  forced  diffusion  is  absent, 

and  when  it  is  further  assumed  that  the  diffusion  is  two-  component  equimolal 

1 1 .. 

counterdiffusion,  then  the  diffusion  equation  simplifies  to  Fick's  second  law’ 
of  diffusion.  When  the  expression  for  the  momentum  flux  (2.6)  iB  substituted 
into  the  equation  cf  motion  (Eq.  (1.  2)  ),  then,  if  the  bulk  viscosity  is  zero,  one 
obtains  the  usual  Navier -Stokes  equations.  When  the  expression  for  the  energy 
flux  (Eqs.  (2.  17)  tp  (2.  19)  ) is  put  into  the  energy  balance  equation  in  terms  of 
the  temperature  (Eq.  (1.9))  one  obtains  the  heat  flow  equation  in  very  general 
form.  When  the  assumption  is  made  that  viscous  dissipation  effects  and  p v - 
effects  are  negligible,  then  the  temperature  equation  for  a pure  fluid  assumes 
a form  which  is  mathematically  analogous  to  Fick's  second  law  of  diffusion. 
Hence  we  see  that  the  equations  of  change  and  the  expressions  for  the  fluxes 
provide  starting  points  for  heat,  momentum,  and  mass  transfer  studies  which 
are  more  general  than  the  usual  starting  equations.  Since  these  equations  and 
expressions  have  all  been  written  in  vector  and  tensor  form,  they  may  be  used 
to  get  the  correct  starting  point  for  calculations  in  cartesian,  cylindrical,  or 


spherical  coprdinates. 
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3.  The  Tranepo rt  Coe f fie  leni £ in  Terms  of  Inter  molecular  Forces 

In  order  to  use  any  results  based  on  solutions  of  the  equations  of  change, 
one  must  know  the  values  of  the  transport  coefficients  involved,  and  in  some 
problems  the  tempi*  r aiure  density,  and  ronrrnfration  dependence  of  these 
quantities  is  also  needed  Whenever  possible,  of  course,  experimental  data 
Should  be  used  it  they  are  available.  Otherwise,  one  should  use  the  best  meth- 
od ot  calculation  which  is  consistent  with  the  accuracy  desired.  In  general,  the 
molt  reliable  methods  are  those  baaed  on  statistical  mechanical  formulae,  which 
relate  the  buik  properties  in  terms  of  the  inter  molecular  forces.  Kence  this  sec- 
tion begins  with  a summary  of  the  various  types  of  forces  between  molecule?,, 
a . Inter  molecular  Force  s and  Potential  Energy  Functions 
In  principle  it  is  possible  to  determine  the  force  between  a pair  of  molecules 
from  a priori  quantum  mechanical  calculations.,  inasmuch  as  the  description  of 
tbe  mechanics  of  any  collection  of  nuclei  and  electrons  la  given  by  the  Schro- 
dinger  wave  equation.  In  practice,  however,  it  is  possible  to  obtain  solutions 
to  the  Schrodinger  equation  for  only  rather  simple  molecular  systems,  and  then 
only  after  making  numerous  approximations  the  physical  reality  of  which  is  dif- 
ficult to  assess.  The  quantum  mechanical  approach  has  nevertheless  provided 
a considerable  amount  of  information  of  a semi -quantitative  character,  which 
has  been  of  great  valu“  in  establishing  tbe  approximate  functional  form  of  the 
interaction  energies. 


An  elementary  discussion  ot  inter  molecular  forces  and  empirical  potential 
energy  functions  may  be  found  in  MTGL  §1-  3,  and  a complete  treatment  is 
given  in  MTGL,  Chs.  12,  13,  and  14.  Sec  also  the  earlier  review  articles 
by  London  (26)  and  Margcnau  (27). 
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In  the  kinetic  theory  formulae  for  the  transport  coefficients  it  ie  the  po- 
tential energy  of  interaction  which  appears  rather  than  the  force  of 

interaction  F(r)  between  the  molecules  in  the  fluid.  For  two  spherical  non- 
polar molecules  the  force  and  potential  energy  are  simply  related  according 
to  the  equation: 

,<*  (30 

F(' ) - - dcp/dr  Tfr)  - Jr  F(r)dr 

Fpr  the  interaction  between  more  complex  mol  *cules,  ihe  potential  energy  io 
a function  not  only  of  the  inter  molecular  separation,  but  also  of  the  angles 
describing  the  mutual  orientation  of  the  molecules.  From  such  an  angle- 
dependent  potential  energy  function  one  can  derive  both  the  force  ar.d  torque 
of  interaction  as  a function  of  tl  a distance  and  orientation. 

It  is  customary,  although  somewhat  arbitrary,  to  divide  intormolecular 
forces  into  two  types  - -short- range  forces  and  long-range  forces.  The  short- 
range  forces  arise  when  the  molecules  come  close  enough  together  for  their 
electron  clouds  to  overlap  appreciably.  These  forces  are  repulsive  and  often 
highly  directional.  There  are  some  experimental  indications  of  the  nature  of 
short  range  inter molec  uiar  torcea  from  crystal  structure,  properties  of  sur- 
faces, and  other  properties,  but  most  of  our  information  comes  from  the  spe- 
cific quantum  mechanical  calculations  which  have  been  made  for  specific  mole- 
cular interactions.  These  calculations  usually  are  not  highly  accurate,  inas- 
much as  many  of  the  integrals  encountered  are  so  complicated  that  they  are 
either  approximated  or  simply  neglected. 
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The  long-range  forces,  on  the  other  hand,  may  he  treated  in  a (airly  rigor- 
our  manner,  f urthermore,  it  is  easy  to  give  formulae  for  various  types  of 
interaction  which  are  applicable  to  a variety  o£  types  o£  molecules.  The  long- 
range  forces  may  be  subdivided  into  electrostatic,  induction,  and  dispersion 
forces.  The  electrostatic  forces  are  those  between  the  various  multipoles  of 
the  molecules  and  are  described  by  electrostatics.  The  induction  forces  are 


those  reculting  from  the  fact  that  a multipole  In  one  molecule  can  induce  a 
multipole  in  the  other  molecule,  which  in  turn  reacts  with  the  multipole  of  the 
first  molecule.  The  dispersion  forces  result,  in  a sense,  from  the  interaction 
of  the  instantaneous  multipoles  which  exist  even  in  neutral  symmetrical  mole- 
cules because  of  the  electronic  motion;  these  forces  cannot  be  described  cor- 
rectly in  terms  of  a simple  phyeical  picture  and  are  derived  by  quantum  me- 
chanical second-order  perturbation  theory. 

For  nonpolar  molecules  the  most  important  long-range  contribution  to  the 
potential  energy  of  interaction  is  the  dispersion  fc  ;e  proportional  to  the  in- 
verse sixth  power  of  the  intermolecular  separation.  The  short-range  force. 


which  is  known  1 be  approximately  exponential  in  form, 


proximate!1/  exponential  in  (enw,  can  be  further  approximated  by  an  inverse 
power  dependence,  usually  r_^^.  It  is  these  considerations  which  form  the 


justification  for  the  Lennar  i-  Jones  (6-12)  potential: 

c$(r)  = 4e  [(cr/r)'2--  (cr/r)fe] 


In  which  <T  is  that  value  of  the  separation  for  which  the  potential  energy  is 
zero,  and  6 Is  the  maximum  attractive  energy.  This  potential  energy 
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I unction  has  been  used  with  considerable  success  for  the  calculation  of  many 
equilibrium  and  non-equilibrium  properties  of  dilute  and  moderately  dense 

£ 

gases  composed  of  molecules  which  are  nonpolar  and  reasonable  spherical. 

Empirical  potential  energy  functions,  such  as  the  Lermard-Jones  potential, 
have  played  a very  *mporfant  role  in  connection  with  the  determination  of 
intermolecular  forces.  Many  physical  properties  can  be  expressed  by  means 
of  statistical  mechanics  in  terms  of  integrals  involving  the  intermolecular 
potential  energy  function.  W lie n an  empirical  potential  function  is  used,  these 
integrals  can  be  evaluated  and  the  results  expressed  in  terms  of  the  parameters 
(such  as  or  and  6 in  the  L.enn«vrd~  Jones  potential)'.  Experimental  data  for  Borne 
bulk  property  can  then  be  used  in  conjunction  with  these  results  io  detei  mine 
the  potential  parameters.  Once  the  potential  parameters  are  known,  other 
physical  properties  can  then  be  calculated.  This  interrelation  of  bulk  proper- 
ties and  intermolecular  forces  is  further  discussed  in  connection  with  Eqs.(3.  7), 
(3.8),  and  (3.9).  Much  more  quantitative  information  about  intermolecular 
forces  has  been  obtained  by  means  of  thiB  semi-empirical  approach  than  by 
direct  quantum  mechanical  calculation. 


* The  following  Lennard -Jones  potential  calculations  are  described  in  MTGL: 
second  and  third  virial  coefficients,  §2.6;  Joule -Thomson  coefficient,  §3  6, 
the  Leonard-Jcnes— Devonshire  equation  of  state  for  dense  gases  and  liquids, 
§4.7;  quantum  calculations  of  the  second  virial  coefficient,  §§6.4c  and  6.  Sc; 
transport  properties  of  dilute  gases,  §8.4  Complete  tables  of  all  tabulated 
functions  are  given. 
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Extensive  calculations  have  recently  been  mads^  for  another  empirical 
function  for  non-polar  moleculeSj  namely  the  modified  Buckingham  (6-exp) 
potential  energy  function.  This  potential  contain*  the  r"^  attractive  term, 
but  the  repulsive  component  is  represented  by  an  exponential  term  thus; 


Cf(r) 


exp  [« (l  - £ 1]  - fe)'6 


(3-1) 


in  -which  there  are  three  adjustable  parameters:  the  maximum  energy  of 
attraction  € . the  separation  corresponding  to  that  energy  rm  , and  the 
parameter  QL  which  is  a measure  of  the  steepness  of  the  energy  of  repul- 
sion.^ This  intermolecular  potential  function  represents  an  improvement 
over  the  Lennard-Jones  ootential  and  has  been  found  to  give  somewhat  better 
agreement  between  calculated  and  experimental  transport  coefficients  and 
virial  coefficients,  particularly  for  and  He. 


t Calculations  of  Mason  (28)  of  the  transport  coefficients  for  the  modified 
Buckingham  potnetial  are  deecribed  in  MTGL,  §8.4.  The  corresponding 
calculations  of  Rice  and  Hirschfelder  (34)  are  described  in  MTGL.,  §3.  7. 
Further  discussion  of  the  modified  Buckingham  potential  has  been  given  by 
Mason  and  Rice  (29). 

' Actually  for  very  small  values  of  the  inter  mole  culaf  separation  x , the 
function  given  in  Eq.  (3.  3)  has  a maximum  value  and  approaches  -minus 
infinity  as  r approaches  aero.  This  would  correspond  to  very  strong 
attraction  at  very  small  separations.  Hence,  when  Eq.  (3.3)  ’s  used  in 
statistical  mechanical  calculations,  one  can  set  CPfr)  * ao  for  values  of 

kM  IK.  - , ' 

rAatLhe  spurious  maximum. 
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b.  The  Transport  Coelficients  of  Dilute  Gas? b^1 

The  earliest  attempts  to  relate  the  transport  coefficients  to  inter  molecular 
forces  were  based  on  the  assumption,  that  molecules  could  be  represented  as 
rigid  spheres  with  diameter  0~^  . Tlvese  early  theories  were  baaed  o«  the 
concept  of  the  "mean  free  path!i--the  average  distance  traversed  by  the  mole- 
cules between  collisions.  A simple  calculation  shows  that  this  quantity  is  giv- 
en  by  L ~ (UT/v2  pjTtr^*).  The  quantity  TTU^fc  appearing  in  the  denominator 
is  called  the  •'collision  cross-section"  for  the  rigid-spherical  molecule.  This 
cross-section  appears  in  all  the  expressions  for  the  transport  coefficient*  and 
represents  the  cross -sectional  area  of  the  imaginary  sphere  surrounding  a mole- 
cule, into  which  the  Renter  of  an  approaching  molecule  cannot  penetrate. 

According  to  the  simple  kinetic  theory  the  transport  coefficients  can  all  be 
expressed  in  terms  of  the  mean  free  path  and  the  average  molecular  velocity 

i ♦ 

v = AffidT/TTm  . Thus  the  coefficients  o.‘  self -diffusion,  viscosity,  and  ther- 

(34) 

(3.5) 

(3.0 


1 Transport  properties  of  dilute  gases  and  binary  gas  mixtures  have  been 

treated  in  a highly  mathematical  form  in  the  book  of  Chapman  and  Cowling(6), 
An  elementary  discussion  of  kinetic  theory  is  given  in  MTGL>,  $1.2,  and 
the  kinetic  theory  of  multicomponent  gas  mixtures  is  given  in  MTGIn  Chs. 

7 and  •?.  Formulae  in  practical  unite  are  given,  and  CO  uupivto  tables  for 
une  in  conjuctiou  with  these  formulae  may  be  found  in  the  appendixes.  A 
summary  of  the  piesenc  status  of  calculating  transport  and  other  proper  - 
ties  of  g?ses  and  liquids  has  been  given  in  a survey  article  written  by  the 
authors  (1). 

^ ~lh ft  d!Scu4*ior  i is  confined  -+c>  Hie  conwdtmium  of  +tv,  proptrfvi  cf  f>ur*  .subiianfi*.  rkrvf  Hy.  irjnitod 
Orfffiow-i  CKvactcrl«i  . » lit  fl(  lef  diffu„Dn,  Wh.s-h  ■>  the.  ,r,<crd. ffiel »n 

OV  psrVlu  of  **k  JO JW  ^ lire.  AHhc*«p  H a Sonrx-.vMt  trhfc.i  .tlinot  who"* 

upp.uMon.  The  mkrd.f  of  ht*vj  isotope;  and  Hu  oliff,i,0n  of  or>Hv,  and  pare 
Art  r>her.or-n».n*  '-huh  fan  b<  de  Py  .'*lf  d'ftvjion  fwrr.iHe . 


mal  conductivity  of  a pure  gas  are  given  by: 

On.  /a  ! irr  I f'  J IT  ro^CT  1 


& = (i/3) 


= (1/3)  f*L  = r 


TTmk'T 


'X  = (1/o)pCvL  * Cx  s/^rnk'!-  ^ 

‘ Tfcj;^ 
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In  which  , and  are  constants;  and  Cv  i»  the  specific  heat 

per  unit  buss,  which  is  (3k/2m)  for  monatomic  substances.  According 
to  the  simple  kinetic  theory  treatment  all  these  quantities  are  (2/3)  . Ac- 
cording to  the  detailed  theory  of  Chapman  and  Cowling  each  of  these  quantities 
is  different:  Cj,  = (3/8)  , = (5/16)  , and  « (25/32)  . This  is  an 

indication  that  the  transport  of  mass,  momentum, . and  energy  cannot  be 
treated  exactly  analogously  to  one  another  as  is  done  in  the  simple  treat- 
ments. Since  for  an  ideal  gas,  j>  = mp/kT  , Eq.  (3.4)  indicates  that  the 
coefficient  of  self -diffusion  should  vary  inversely  to  the  pressure  and 
directly  as  the  3/2  power  of  the  temperature.  According  to  EqB.  (3.5)  and  (3.6) 
the  coefficients  of  viscosity  and  thermal  conductivity  of  a dilute  gas  should 
be  independent  of  the  pressure  and  directly  proportional  to  the  square  root 
of  the  temperature.  Thus  the  simple  theory  predicts  the  correct  pressure 
dependence  of  the  dilute  gas  transport  coefficients  and  comes  close  to  giving 
the  correct  temperature  dependence.  It  also  provided  the  early  kinetic 
theory  workers  with  a means  for  estimating  moleculav  diameters  from 
bul*.  properties. 

The  rigorous  kinetic  theory  developed  by  Chapman  and  Ennkog  and 
described  in  the  monograph  of  Chapman  and  Cowling  (6)  gives  an  indication 
as  to  how  the  results  of  the  simple  kinetic  theory  havo  to  he  modified  in 
order  to  apply  to  real  molecule*  rather  than  rigid  spheres.  The  rigorous 
theory  does  not  An  the  course  of  its  development  employ  ths  concept  of  the 
mean  free  path.  The  rigorous  theorv  has  us  its  starting  point  the  Boltswnann 
integro  -differential  equation  for  the  velocity  distribution  function.  Although 
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the  details  of  the  matuematlcal  development  are  quite  intricate  the  results 
for  pure  substances  may  be  expressed  in  -a  simple  form  reminiscent  of 
the  result*  given  above  for  the  simple  kinetic  theory.  Thus  the  rigorous 
expressions  for  the  first  approximations  to  the  coefficients  of  self-diffusion, 
viscosity,  and  thermal  conductivity  for  a dilute  monatomic  gas  are: 

[»],  ' 

/ 1 

(3.©) 


A.  ^1Trn^  -! 

6 f 

tri,  - I 


Kb  ircr*^ 


[>1  - =5l  JJU^L  r 

LAJl  u 


J5  k. 

4- 


few, 


(3.9) 


in  which  the  Q,  -functions  indicate  the  deviation  of  real  gas  L«*havior  from 
rigid  Bphere  behavior. f 

The  cCI  -functions  are  complicated  integrals  which  depend  upon  the 
force  law  between  the  molecules.  When  the  L,ennard-Joneo  potential  is 
used,  cr  in  Eqs.  (1,7),  (3.8),  and  (3.9)  is  the  a ir-  the  Lennard-Jones 
potential.  The  -functions  are,  then,  functions  of  a.  dimensionless 
temperature  kT J6  , in  which  the  <■  is  the  maximum  energy  of  attrac- 


In  terras  of  the  notation  used  in  MTGL: 

Actually  there  is  a whole  set  of  in  terms  of  which  the  transport 

properties  can  be  expressed.  The  ones  with  indices  (1,1)  and  (2,  2)  are 
the  most  important,  hut  additional  -function*  are  required  for  mixtures 
and  for  higher  approximations,  All  of  the  -functions  needed  for  calcula- 
tions may  be  found  in  MTGL.  (18)  for  the  following  potential  energy  functions: 
rigid  spheres,  point  centers  of  repulsion,  the  Sutherland  model,  the  square 
well,  the  .Lennard-Jones  (6-12)  potential,  and  a modified  Buckingham  (£>-exp) 
potential. 
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tion  appearing  in  the  Leonard -Jones  potential  function.  Hence  the  transport 
coefficients  are  expressed  as  universal  functions  of  a reduced  temperature 
and  depend  parametrically  on  the  constants  CT  arid.  £.  which  are  characteristic 
of  the  substance  under  consideration.  Hence  if  a*  and  e are  known  for  some 
substance,  the  self -diffusion,  viscosity,  or  thermal  conductivity  can  be  calcu- 
lated from  the  above  equation  by  use  of  the  tabulations  available  for  the  if! 
functions.  Or,  if  the  coefficient  of  viscosity,  for  example*,  is  known  -for  emit 
substance  at  two  different  temperatures,  then  or  and  e for  that  substance  can 
easily  be  found  by  the  solution  of  two  simultaneous  equations,  This  illustrates 
the  way  in  which  inter  molecular  force  information  is  obtained  from  the  measur- 
merits  of  bulk  properties. 

Eqs.  (3.7.  8;  9)  ar*  based  cn  the  rigorous  kinetic  theory  for  monatomic 
gases.  The  use  of  F.qs.  (3.  7,  8)  for  diffusion  and  viscosity  of  paly ato mi?,  gases 
turns  out  to  be  quite  a good  approximation,  provided  that  the  molecules  are  not 
polar  and  t&ot  too  elongated.  Eq.  (3.9)  for  thermal  conductivity,  however,  must 
be  modified  for  polyatomic  molecules  by  means  of  the  approximate  “Eucken  cor- 
rection" thus: 


Until  the  further  development  of  the  existing  theories  for  the  transport  phenomena 
of  polyatomic  molecules  (42,  10).  the  Euc ken  Correction  can  be  used  to  a 
reasonably  good  approximation,' 


| A theoretical  justification  tor  tjtt  v.uckea  correction  is  given  in  MTGL,  §7.  6b. 


WJfi-ONR-7 
i«  May  1954 
7 -3'- 

Besidea  formulae  and  tables  for  computing  the  first  approximation  to  the 
transport  coefficients  of  pure  substance*  as  discussed  above,  there  are  also 
tabulated  functions  available  for  the  Lennard-Jensn  potential  fat  us*  with  formu- 
lae for  the  transport  coefficients  of  binary  and  multicomponent  aal?*turee  and  for 
the  higher  corrections  ts  the  transport  coefficients.  The  formula*  are  <quite 
lengthy  and  are  not  given  here  since  they  are  given  in  full  elsewhere.  Generally 

speaking  these  results  are  easy  to  use  and  give  excellent  agreement  with  the  ex- 

(U8) 

pcrimental  data.  Quantum  effects  have  boon  shown  to  be  negligible  at  room 
temperature  and  above,  even  for  helium  and  hydrogen  IS). 

The  situation  with  respect  to  dilute  gases  is  then  quite  satisfactory  at  the 
present  time,  particularly  as  far  as  nonpolar  molecules  are  concerned.  Much 
remains  to  be  done  on  Ihc  development  of  the  kinetic  thepnry  of  polyatomic  g^.ses 
and  particularly  the  applications  to  polar  gases.  Aside  from  the  very  general 
and  highly  mathematical  developments  of  Wfcng  Chang  and  Uhlenbeck  (42)  and 
de  Boer  (10),  the  only  other  work  which  has  been  done  on  the  kinetic  theory  of 
non-spherical  molecules  has  been  that  on  several  special  models,  'such  as  per- 
fectly rough  spheres  (5,  32),  loaded  spheres  (23),  and  rigid  ellipsoids  (22). 

c.  Transport  Coefficients  of  Dense  Gases  and  Llqujdst 

At  the  present  time  there  are  f^ur  main  approaches  used  in  the  study  of  the 
transport  properties  pi  dense  gases  and  liquids.  None  of  these  methods  can  be 
regarded  as  entirely  satisfactory,  and  indeed  none  of  them  rc-slly  provides  usa- 
ble results  as  far  as  interrelating  bulk  properties  with  intermoleculsr  forces 
is  concerned.  describe  these  methods  in  order  of  increasing  mathematical 

CO  alp  Ic  xit  y . 
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Th»  principle  af  cor  responding  states,  which  has  bean  of  great  value  for 
estimating  ^p-V-T  and  thermodynamic  properties  of  fluids  (20),  can  al>o 
be  applied  to  the  correlation  of  transport  coefficient  data.  That  is,  the-  tram*- 

4- 

port  cssffients  divided  by  their  values  at  the  critical  point,  mr.y  be  expressed 
as  universal  functions  of  the  reduced  pressure  and  reduced  temperature: 


“ 

- &r.(pr,Tf) 

(3,11) 

|^r  ** 

^/Mc 

« \ir(  prX) 

(3.12) 

Ar  *» 

a / 7\c 

“ Mr  tpr.^r) 

(3.13) 

in  which  * T/T„  and  pr  = p/pe.  Equation  (3.12)  provides  the  basis 
for  Uyehara  and  Watsbn's  (41,  20)  generalized  viscosity  chart,  and  Eq.  (3.  13) 
for  Gammon's  (14)  generalized  thermal  conductivity  chart.  The  latter  was 
prepared  with  relatively  few  experimental  data.  Aa  soon,  as  more  data  on  high 
density  transport  coefficients  are  available,  they  should  be  correlated  in  the 
form  of  generalized  charts  based  on  the  principle  of  corresponding  states. 


^ The  material  presented  in  this  section  is  treated  more  extensively  in  LvfGL, 
CH.  9. 

4 Other  methods  of  reducing  the  variables  have  been,  suggested.  One  can,  for 
example,  reduce  the  transport  coefficients  by  dividing  them  by  the  approp- 
riate combination  of  critical  constants,  thus: 

H"  V kTt/  n3^4' fA  » fA*(fy,T-)  (3,12.*) 

la  etmihs?  mathod  one  divides  the  transport  coefficient  by  lit  limiting  aero- 
prea:*ure  vzlue,  thus; 

r-  - M.#(pr>Tr)  (3.12b) 

This  seethed  was  used  by  Comings,  Mayland,  and  Egly  (7)  in  the  preparation 
of  a generalised  chart. 
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*31  The  theory  of  rate  processes  ac  developed  by  Eyring  and  hi«  collaborators  (15) 
has  been  moderately  successful  in  explaining  the  transport  phenomena  in  liquids 
on  a simple  pictorial  baeie.  Such  an  approach  does  not  lead  to  expressions  for 
the  transport  coefficients  in  terms  of  inter  molecular  forces,  but  rather  to  re- 
lations between  various  macroscopic  quantities; 


a a 

(b  = nb  e*p  (0.406  ^~%/RT) 

(3.14) 

= nh  exp  (3,8  Tb/T ) 

y 

(3. 14a) 

£>  = n'/3kT/K- 

(3.15) 

A - 787?  k r,ft  c 

T rm/wbmic 
[ liquids 

(3  1 fe) 

7\  = 2.80  k n“'J  X~'k  c 

j pol^atomsc 
L iio|Liidi 

(3.»7) 

In  these  relations  h is  Planck's  (constant  and  c la  the  speed  of  sound  in  the 
liquid.  Eq.  (3.  14a)  is  a simplified  form  of  Eq.  (3.  14)  made  possible  by  the 
use  of  Trouton's  rule.  Equation  (3.  16a)  is  a generalisation  of  Eq.  (3.  16) 
made  possible  by  the  use  of  a modified  Eucken  correction.  The  Eyring  theory 
has  been  used  to  study  the  transport  properties  of  liquida  composed  of  nonspheri- 
cal  molecules  such  as  long  chain  hydrocarbons  and  high  polymerc  , It  ie  the  only 
theory  of  transport  phenomena  whie*>  predicts  non -newt?  ui  an  flow. 

The  kinetic  theory  for  rigid  spheres  developed  by  Enakog,  ( 6 ^ 18)  is  probably 
the  best  theory  available  at  the  present  time  fsr  describing  the  transport  coef- 
ficients of  dense  gases.  The  big  difficulty  in  the  development  of  a kinetic  theory 
for  condensed  systems  is  that  one  must  understand  certain  aspect#  of  three  - 
molecule  and  higher  order  collisions.  For  the  rigid  sphere  model,  however, 
it  is  theoretically  impossible  for  three  molecules  tc  collide  at  exactly  the  same 


WJ&-ONR-7 
12  May  1954 

-40- 


meraent.  Hence  a complete  theory  may  be  developed  by  considering  two-body 
collisions  only.  The  final  results  of  RnSkeg  theory  may  be  Summarised  in 
terms  of  the  following  relations,  which  interrelate  the  reduced  transport  coef- 
ficients (based  on  sere  pressure  values)  and  the  compressibility  factor  for 


rigid  sphsrei: 

(P'/^XV/bo)  - (1/y)  4 0.8  + 0.7fcly 

(3.17) 

(K/R "*)(ty/b*)  « 1.0 QZ  y 

0 18) 

(WA°)Ctf/bo)  = (i/y)  + i-Z  -b  0.753  y 

(3.19) 

(£>/4b°)(V/b o)  « (l/y) 

(3.20)  f 

In  these  relations  bQ  * (2/3)ttNor^3  and 

j « (bV/RT)  - 1 - (bo/^)  + 0-62.50  Cb^/W^ 

fo  P;N. 

+ a 286?  (W/Vf  + 0.115  (b„/Vf+  - 


which  ia  obtained  from  the  virial  equation  of  state  for  rigid  spheres.  Although 
these  results  were  obt*i.»ed  for  rigid  spheres,  Enskog  showed  that  for  the  availa- 
ble experimental  data  the  following  procedure  could  be  used,  fot  teal  gases;  One 
replaces  the  pressure  p i \ Kq.  (3.  21)  by  the  "theru&l  pressure*'  T (dp/3T),^. 
co  that  y may  be  determined  from  the  experimental  p-'V-T  data  from  the 
rolationt 

J - (WR  A^P/^V  ~ i (3,22) 

( For  rigid  spheres  Eqs.  (3.  21)  and  (3.  21a)  are  the  fcaic'4,  but  such  is  not  the 
ctse  for  real  molecules. ) Then  one  fits  the  minimum  in  iiv*  curve  of  the  vis- 
cc-oity  versus  y in  order  to  specify  a value  of  bQ  . It  will  be  interesting  to 

^ may  be  calculated  as  oC  = MD where  oo1  is  the  coefficient  of  eeif- 

di/iusion  calculated  at  1 aim.  pressure  according  to  E?,.  (4.7). 
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teat  this  method  further  when  more  experimental  data  become  available . The 

Enskog  iluvory  has  been  extended  Lo  binary  mixtures  by  Thorne  4<c)* 

The  -methods  of  non- equilibrium  statistical  mechanic*  have  been  applied  to 

(21; 

the  problem  of  transport  phenomena  in  dense  systems  by  Kirk  woodland  bv  Born 

r \ 

and  Green  <3),  The  firial  results  are  given  in  terms  of  n no>ne%uJUibriuaa  radial 
distribution  function.  This  function  is  the  solution  of  an  equation,  whtch, 

however,  involve#  the  next  higher  order  distribution  function.  A M“ ups r position 
approximation"  is  introduced  in  ordv»r  to  get  rid  of  the  higher-order  distribution 
function,  but  the  validity  c£  this  assumption  has  not  been,  fully  assessed.  At  the 
present  time  only  limited  calculations  have  been  made,  and.  hence  the  method 
do«£  not  yet  provide  a means  for  practi  cal  computation.  It  is  hoped  that  this 

approach  will  ultimately  lead  to  the  calculation  of  dense  gns  properties  with 

\ V' 

accuracy  comparable  to  that  for  dilute  gases  according  to  the  Ghapman-Enskog 
theory. 


This  work  was  carried  out  under  Contract  K7o:*?-28511,  Office  of  Naval 
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Notation 


Rivman  Letters 


3^  * activity  of  species  ,1i,c 

x phenomenological 
coefficients 

tb  * (2/2)JTNctr* 

Qy  * specific  heat  at  con- 
stant volume 

dt  x driving  force  defined 
J in  Ep.  1.  10 


T\/r\t- 


•W  ±J\i 


s substantial  deri  vat)  ve 
operator 


multicomponent  dif- 
fusion coefficients 


x.  binary  diffusion  coef- 
* ficients 

D?  x multicomponent  thermal 
diffusion  coefficients 


= external  force  acting  on 
species  i 


3 

6 


x rate n£  irre.versible  pro- 
duction o£  enltropy 
e Gibbs  free  energy 


H 

i 

ji 


» enthalpy 

x unit  tensor 
x ma«*  flux  of  species  i 
' (gm/cm2  sec)  with  re- 
spect to  mass  average 
velocity,  v 


Js 


x mass  flux  of  secies  i 
fmoies/cm2  sec)  'cith 
respect  to  molar  aver- 
age velocity 


Jh  = fl voces  in  the  thermo- 
dynamics of  irrovers- 
ibi*  thermodynamic  s 

k 3 Boltsmana'e  constant 

kT  = thermal  diffusion  ratio 

Kl  “ rate  of  production  of 

species  f by  chemi- 
c si  xs  action  (gm/sec) 

L x mean  free  path  in  kinetic 
theory 

rrii  x mol«cul4x  DttSOS  of  species  i 

n x number  density  {number  of 

molecules  per  unit  volume) 

Hi  = number  density  of  species  i 

nt  = number  of  chemical  species  in 
multicomponent  system 

nr  = numLer  of  reactions  in  multi- 
eetwpSitgnt  system 

Ni  = mass  flax  of  i (molee/cm2  sec) 
stith-TO-spacJ  to  stationary  axes 

FJ  x A vo gad vo's  number 

p = static  pressure 

= partial  pressure  of  specie#  i 

p = pressure  te?*c'*r 
yjS* 

P » vinous  pari  pi  pressure  tenoor 

C|  - energy  ftp*  v»etor 
c(W  _ beat  fiiM  vedrt 


WIS-ONR-7 
12  May  1954 
Notation 


r * 

distance  between  two  col- 

k. 

=■  coefficient  of  bulk  vistfOpity 

liding  molecules 

> 

=•  coefficient  of  thermal  con- 

G  = 

rate  of  the  sth  chemical 

ductivity 

reaction 

Aa 

f <■ 

R . 

= driving  force  defined  in  Eq. 

gas  constant 

1.  13 

S , 

entropy 

H- 

= coefficient  of  viscosity 

t e 

time 

\k 

- chemical  potential  of  i1*1 

-r 

1 * 

absolute  temperature 

i 

species  (permolacule) 

\ • 

TT 

= 3. i4i59... 

normal  boiling  tempera- 

ture 

P 

= mass  density  of  fluid 

r t 

o = 

internal  energy 

Pi 

= mass  density  of  il“  rom- 

pc  no  tit  in  fluid 

•v  = 

mas*  average  velocity 

• . t i 

(stream  velocity) 

cr 

= parameter  in  inter  mo  lecular 

V - 

potential  energy  function 

average  molecular  speed 

°R 

- diameter  of  rigid  sphere 

V . 

volume 

molecule 

Yld>« 

diffusion  velocity 

cr 

- entropy  flux 

Xi  = 

mol  fraction  of  species  i 

c«r) 

= intermolecular  potential 

h*~  coordinate --X.  v.  or  z 

energy  function 

*k.  * 

iiii  ] 

1 

xh  - 

,r 

y * functions  of  the  reducer 

generalised  forces  in  the 

« j 

n»*  i 

1 temperature  which  are 

thermodynamics  of  irre- 
versible thermodynamics 

u nX  ■' 

introduced  in  Eqs.  3.  7,  8,  9 

Subscripts 

y a 

quantity  defined  in  §3c 

a chemical  species  present  in 

M>k 

X * 

chemical  affinity  for  the  . 

mvdticOHJpOnent  mixture 

*th  reaction 

5 

- quantities  pertaining  to  sth 

Greek 

letters 

c he  mic  al  re  action. 

«6i  * 

*re si Stance  factors" 

C 

* critical  quantities 

€ a maximum  energy  of  attrac- 
tion in  into r molecular 
potential  energy  function 


r - qu&htitiea  reduced  by  division 
by  the  corresponding  critical 
quantities 
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Superscripts 

O *-.  ideal  gas  properties 

"I"  * transpose  of  a tensor  or 

dyadic  product 

Above  Symbols 

A = quantity  per  unit  ma3s 

r-*/  = quantity  per  mole 

— = partial  molal  quantity 
Belov/  Symbols 

— = vector  quantity 
= tensor  quantity 

Vector  and  Tensor  Operations 

1,  * unit  vectors  in  x,  y,  x, 

- directions 

2?  = tVn  + jVy  + 

VT  » L(bT/bx)  + im/by)  + k (bT/ta) 

(V*a)  = /bx)  -t-  (bx^/by)  + (bir«  /££) 

(Vp)x*  (bp^/ax)  + (bp^/dy)  + (bp^/br) 

(u-V)T  s (bT/bx)  4 u^(bT/by)  4 v*n(b7Ybz) 

(tJtVa)  * K ('a’S,/dr)  + P^(b\r,/ty)  + tv,  (bir»/eo.) 

+ ^(bi^/by)  4 p^(buy/c)z) 

+ hx(^r/^<)  + ^(buv/by)  i ^(b^/bz) 

DT/Dt  . "bT/bt  4 (V*V)T 


£?£.FEfcffN£?S 
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